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Abstract — Despite the central role of mobility in wireless 
networks, analytical study on its impact on network performance 
is notoriously difficult. This paper aims to address this gap by 
proposing a random waypoint (RWP) mobility model defined 
on the entire plane and applying it to analyze two key cellular 
network parameters: handover rate and sojourn time. We first 
analyze the stochastic properties of the proposed model and 
compare it to two other models: the classical RWP mobility model 
and a synthetic truncated Levy walk model which is constructed 
from real mobility trajectories. The comparison shows that 
the proposed RWP mobility model is more appropriate for 
the mobility simulation in emerging cellular networks, which 
have ever-smaller cells. Then we apply the proposed model 
to cellular networks under both deterministic (hexagonal) and 
random (Poisson) base station (BS) models. We present analytic 
expressions for both handover rate and sojourn time, which have 
the expected property that the handover rate is proportional 
to the square root of BS density. Compared to an actual BS 
distribution, we find that the Poisson- Voronoi model is about as 
accurate in terms of mobility evaluation as hexagonal model, 
though being more pessimistic in that it predicts a higher 
handover rate and lower sojourn time. 

I. Introduction 

The support of mobility is a fundamental aspect of wireless 
networks JT), G) ■ Mobility management is taking on new 
importance and complexity in emerging cellular networks, 
which have ever-smaller and more irregular cells. 

A. Background and Proposed Approach 

As far as cellular networks are concerned, typical questions 
of interest include how mobility affects handover rate and 
sojourn time. Handover rate is defined as the expected number 
of handovers per unit time. It is directly related to the network 
signaling overhead. Clearly, the handover rate is low for large 
cells and/or low mobility, but smaller cells are necessary to 
increase the capacity of cellular networks through increased 
spectral and spatial reuse. Thus, analytic results on handover 
rate will be useful for network dimensioning, and in order 
to understand tradeoffs between optimum cell associations 
and the undesired overhead in session setup and tear downs. 
Sojourn time is defined as the time that a mobile resides 
in a typical cell. It represents the time that a BS would 
provide service to the mobile user. In the case of a short 
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sojourn time, it may be preferable from a system-level view 
to temporarily tolerate a suboptimal BS association versus 
initiating a handover into and out of this cell. This concept 
is supported in 4G standards with a threshold rule [3 j, but it 
is fair to say that the theory behind such tradeoffs is not well 
developed. Though handover rate is inversely proportional to 
expected sojourn time, their distributions can be significantly 
different, which motivates us to study them separately in this 
paper. 

To explore the role of mobility in cellular networks, par- 
ticularly handover rate and sojourn time, mobility modeling 
is obviously a necessary first step. In this paper, we focus 
on RWP mobility model originally proposed in [4| due to its 
simplicity in modeling movement patterns of mobile nodes. 
In this model, mobile users move in a finite domain A. At 
each turning point, each user selects the destination point 
(referred to as waypoint) uniformly distributed in A and 
chooses the velocity from a uniform distribution. Then the user 
moves along the line (whose length is called transition length) 
connecting its current waypoint to the newly selected waypoint 
at the chosen velocity. This process repeats at each waypoint. 
Optionally, the user can have a random pause time at each 
waypoint before moving to the next waypoint. In this classical 
RWP mobility model, the stationary spatial node distribution 
tends to concentrate near the center of the finite domain and 
thus may be inconvenient if users locate more or less uniformly 
in the network [5 1. Another inconvenience is that the transition 
lengths in the classical RWP mobility model are of the same 
order as the size of the domain A, which seems to deviate 
significantly from those observed in human walks JSJ. 

To solve the inconveniences mentioned above, we propose 
a RWP mobility model defined on the entire plane. In this 
model, at each waypoint the mobile node chooses 1 ) a random 
direction uniformly distributed on [0, 2ir], 2) a transition length 
from some distribution, and 3) a velocity from some distribu- 
tion. Then the node moves to the next waypoint (determined by 
choice 1 and 2) at the chosen velocity. As in the classical RWP 
mobility model, the node can have a random pause time at 
each waypoint. Note that human movement has very complex 
temporal and spatial correlations and its nature has not been 
fully understood Q, |8). It is fair to say that none of the 
existing mobility models are fully realistic. The motivation of 
this work is not to solve this open problem. Instead, we aim 
to propose a tractable RWP mobility model, which can be 
more appropriate than the classical RWP mobility model for 
the study of mobility in emerging cellular networks, so that it 
can be utilized to analyze the impact of mobility in cellular 
networks to provide insights on network design. 
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After presenting the proposed mobility model, we analyze 
its associated stochastic properties and compare it to two other 
models: the classical RWP mobility model and a synthetic 
truncated Levy walk model J6] which is constructed from real 
mobility trajectories. The comparison shows that the proposed 
model is more appropriate for the mobility simulation in 
emerging cellular networks. Then the proposed RWP mobility 
model is applied to cellular networks whose BSs are modeled 
in two conceptually opposite ways. The first is the traditional 
hexagonal grid, which represents an extreme in terms of 
regularity and uniformity of coverage, and is completely 
deterministic. The second is to model the BS locations as 
drawn from a Poisson point process (PPP), which creates a set 
of BSs with completely independent locations J5), iflOl . Under 
the PPP BS model, the cellular network can be viewed as a 
Poisson- Voronoi tessellation if the mobile users are assumed 
to connect to the nearest BSs. As one would expect, most 
actual deployments of cellular networks lie between these 
two models, both qualitatively - i.e. they are neither perfectly 
regular nor perfectly random - and quantitatively, i.e. the SINR 
statistics and other statistical measures are bounded by these 
two approaches [ 1 1 1. Thus, both models are of interest and we 
apply the proposed RWP mobility model to cellular networks 
under both models. Analytic expressions for handover rate and 
sojourn time are obtained under both models, some of which 
are quite simple and lead to intuitive interpretations. 

In some aspects the proposed RWP mobility model is simi- 
lar to the random direction (RD) mobility model. Though the 
classical RD mobility model does not have a non-homogenous 
spatial distribution as in the classical RWP mobility model, 
its transition lengths are still of the same order as the size 
of the simulation area, which seems to deviate significantly 
from those in human walks \12\. This inconvenience may be 
solved by a modified RD mobility model [13], in which the 
transition length distribution is co-determined by the velocity 
distribution and moving time distribution. This may complicate 
the theoretical analysis in this paper. This motivates us to 
propose a mobility model with transition length distribution 
directly specified. Besides, transition length data seems to be 
more readily available than moving time data in many real 
data sets 0, |8). Thus, it might be easier to use the proposed 
RWP mobility model when fitting the mobility model to real 
mobility trajectories. 

Before ending this subsection, it is worth mentioning 
the many trace-based mobility models lfl4l - ll24l . One main 
drawback of trace-based mobility models is that, due to the 
differences in the trace acquiring methods, sizes of trace 
data, and data filtration techniques, mobility model built on 
one trace data set may not be applicable to other network 
scenarios. Moreover, trace-based mobility models are often 
not mathematically tractable (at least very complicated), pre- 
venting researchers from analytically studying the performance 
of various protocols and/or getting quick informative results in 
mobile networks. In contrast, random synthetic models, includ- 
ing random walk, the classical and the proposed RWP, Gauss- 
Markov, etc., are generic and more mathematically tractable. 
In addition, some aspects of some random synthetic models 
(including the proposed RWP model) can be fitted using traces, 



as done in |[T9l . Nevertheless, trace-based mobility models 
are more realistic and scenario-dependent, which is crucial to 
perform reliable performance evaluation of mobile networks. 
Hence, both random and trace-based synthetic models are 
important for the study of mobile networks. In particular, when 
evaluating a specific protocol in mobile networks, researchers 
can utilize random synthetic models for example the proposed 
RWP model to quickly get informative results. Meanwhile, 
researchers can build specific scenario-dependent trace-based 
model by collecting traces if possible. Then based on the 
constructed trace-based model researchers can further perform 
more extensive and reliable evaluation of the protocol beyond 
just using random synthetic models. 

B. Related work 

The proposed RWP model is based on the one originally 
proposed in [4 |. Due to its simplicity in modeling movement 
patterns of mobile nodes, the classical RWP mobility model 
has been extensively studied in literature 0, Il25l - ll27l . These 
studies analyzed the various stochastic mobility parameters, 
including transition length, transition time, direction switch 
rate, and spatial node distribution. When it comes to applying 
the mobility model to hexagonal cellular networks, simulations 
are often required to study the impact of mobility since the 
analysis is hard to proceed 0, 0. Nonetheless, the effects of 
the classical RWP mobility model to cellular networks have 
been briefly analyzed in [26] and a more detailed study can 
be found in ll28l . However, as remarked above, the classical 
RWP model may not be convenient in some cases. In contrast, 
we analyze and obtain insight about the impact of mobility 
under a hexagonal model through applying the relatively clean 
characterization of the proposed RWP model, as [28 1 did 
through applying the classical RWP model. 

The application of the proposed RWP mobility model to 
cellular networks modeled as Poisson- Voronoi tessellation re- 
quires stochastic geometric tools, which are becoming increas- 
ingly sophisticated and popular ll29l - lf32l . As far as mobility 
is concerned, [33] proposed a framework to study the impact 
of mobility in cellular networks modeled as Poisson- Voronoi 
tessellation. In particular, the authors proposed a Poisson line 
process to model the road system, along which the mobile 
users move. Thus, the mobility pattern in 11331 is of large scale 
while the RWP mobility model in this paper is of small scale. 
Our study can be viewed as complementary to 11331 . For exam- 
ple, our result indicates that if cells decrease in size such that 
the BS density per unit area is increased by 4 times, then the 
handover rate would be doubled. Note that recent work showed 
that the PPP model for BSs was about as accurate in terms of 
SINR distribution as the hexagonal grid for a representative 
urban cellular network ifTTI . Interestingly, we find that this 
observation is also true for mobility evaluation, though the 
Poisson- Voronoi model yields slightly higher handover rate 
and lower sojourn time and thus is a bit more pessimistic than 
the hexagonal model (which is correspondingly optimistic). 

We briefly summarize the contributions of this work: 1) We 
propose a tractable RWP mobility model which overcomes 
some inconveniences of the classical RWP mobility model and 
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is more appropriate for mobility study in cellular networks. 
2) We obtain analytical results for handover rate and sojourn 
time in cellular networks. Though a specific transition length 
distribution is assumed, the analysis in this paper is quite 
general and can be extended to any other transition length 
distribution which has finite mean. 3) We connect the mobility 
results for the two conceptually opposite models of cellular 
networks: hexagonal and Poisson-Voronoi tessellation. 

II. Proposed RWP mobility model on infinite 

PLANES 

We describe the proposed RWP mobility model in this 
section. As in the description of the classical RWP model 
(see, e.g., 0), the movement trace of a node can be 
formally described by an infinite sequence of quadruples: 
{(X n -i, X n , V n , S n )}n£fi , where n denotes the n-th move- 
ment period. During the n-th movement period, X n -i denotes 
the starting waypoint, X n denotes the target waypoint, V n 
denotes the velocity, and S n denotes the pause time at the way- 
point X n . Given the current waypoint X n -\, the next way- 
point X n is chosen such that the included angle between the 
vector X n — X n _i and the abscissa is uniformly distributed 
on [0,27r] and the transition length L n —\\ X n — X„_i || is 
a nonnegative random variable. The selection of waypoints is 
independent and identical for each movement period. 

Though there is a degree of freedom in modeling the 
random transition lengths, we focus on a particular choice in 
this paper. Specifically, the transition lengths {Li,L2, ■■■} are 
chosen to be independent and identically distributed (i.i.d.) 
with cumulative distribution function (cdf) 

P(L < I) = 1 - cxp(-A7rZ 2 ),Z > 0, (1) 

i.e., the transition lengths are Rayleigh distributed. Velocities 
V n are i.i.d. with distribution Py(-). Pause times S n are also 
i.i.d. with distribution Pg(-). This selection bears an interesting 
interpretation: Given the waypoint X n -i, a homogeneous PPP 
<£>(n) with intensity A is independently generated and then the 
nearest point in $(n) is selected as the next waypoint, i.e., 
X n = argmin a , e4 ,(„) || x - X n _i || . 

Under the proposed model, different mobility patterns can 
be captured by choosing different A's. Larger A statistically 
implies that the transition lengths L are shorter. This further 
implies that the movement direction switch rates are higher. 
These mobility parameters may be appropriate for mobile 
users walking and shopping in a city, for example. In contrast, 
smaller A statistically implies that the transition lengths L are 
longer and the corresponding movement direction switch rates 
are lower. These mobility parameters may be appropriate for 
driving users, particularly those on the highways. This intuitive 
result can also be observed in Fig. Q] which shows 4 sample 
traces of the proposed RWP model. 

III. Stochastic properties of the Proposed RWP 

MOBILITY MODEL 

In this section, we first study the various stochastic proper- 
ties of the proposed RWP mobility model. Stochastic prop- 
erties of interests include transition length, transition time, 
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Fig. 1. Sample traces of the proposed RWP mobility models. The transition 
lengths are statistically shorter with larger mobility parameter A, and vice 
versa. 

direction switch rate, and the spatial node distribution. We then 
perform simulation to compare the proposed RWP mobility 
model to the classical RWP mobility model and a synthetic 
Levy Walk model proposed in (6|, which is constructed from 
real mobility trajectories. 

A. Transition length 

We define transition length as the Euclidean distance be- 
tween two successive waypoints. In the proposed model, the 
transition lengths can be described by a stochastic process 
{L ra }„ e N where L n are i.i.d. Rayleigh distributed with 



Note that the transition lengths are not i.i.d. in the classical 
RWP mobility model. Indeed, a node currently located near the 
border of the finite domain tends to have a longer transition 
length while a node located around the center of the finite 
domain statistically has a shorter transition length for the next 
movement period. As a result, it is difficult to obtain the 
probability distribution for each random transition length. 

Nevertheless, the random waypoints X n are i.i.d. in the 
classical RWP mobility model, which is obvious since they 
are selected uniformly from a finite domain and independently 
over movement periods. This property forms the basis for the 
analysis of the classical RWP mobility model (see, e.g., 0, 
[25 1). In contrast, the waypoints in our proposed model are 
not i.i.d. but form a Markov process. 

B. Transition time 

We define transition time as the time a node spends during 
the movement between two successive waypoints. We denote 
by T n the transition time for the movement period n. Then 
T = L/V where we omit the period index n since T„ are i.i.d.. 
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Denote Velas the range of the random velocity V. Given 
any velocity distribution -Py(-), the probability distribution of 
T is given as follows. 

Proposition 1. The cdfof the random transition time T is given 
by 



P(T <t) = l- [ cxp(-\nv 2 t 2 ) dP v (v),t > 0. 
Jv 



(3) 



The proof of Prop. Q] is omitted for brevity. As a specific 
application of Prop. Q] the following corollary gives closed 
form expressions for transition times under two types of 
velocity distributions. 

Corollary 1. 1) If V = v where v is a positive constant, the 
pdf of transition time T is 

f T (t) = 2 7 ^A^/ 2 ^e- Wt, 



2) If V is uniformly distributed on [v 
transition time T is 



,t>0. (4) 

mm, Vmax], the pdf of 
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Next, we derive the mean transition time. Instead of apply- 
ing the cdf of T, we notice that 



E{T\ = J5[|] = E\L\E\^\ 



y obtain 

E[T] = 



From (0, we can easily obtain that, if V = v, 

1 

and, if V is uniformly distributed on [v m in,v n 
In 



E[T] 



(6) 



(7) 



(8) 



2v'A(w max — Wmin) 

From (0, it is clear that v mul > is required to ensure finite 
expected transition time. 

C. Direction switch rate 

In this subsection, we study the direction switch rate, which 
is the inverse of the time between two direction changes and 
thus characterizes the frequency of direction change and is 
also a mobility parameter of interest [5 1. To this end, we first 
introduce the period time, the time a node spends between two 
successive waypoints. In particular, period time T p is the sum 
of the pause time and the transition time, i.e., T p = T + S, 
where S denotes the random pause time. Then the direction 
switch rate denoted by D is defined to be the inverse of 
the period time: D = 1/T p , whose probability distribution 
is described in the following proposition. 

Proposition 2. The cdf of the direction switch rate D is given 



by 

P{D <d) = 



I ^ 2( d - s f 

exp ~\irv 



s Jv 



d 2 



dP v {v) dP s (s), 
(9) 



where 5 £ R denotes the range of the random pause time. 

We omit the proof of Prop. [2] for brevity. Given the distri- 
butions of velocity and transition time, i.e., Pv(v) and Ps(s), 
the distribution of direction switch rate D can be found by 
Prop, m 



D. Spatial node distribution 

In this subsection, we study the spatial node distribution. To 
this end, let Xo and Xi be two successive waypoints. Given 
Xq, we are interested in the probability that the moving node 
resides in some measurable set A during the movement from 
X to X\. We first derive the spatial node distribution given 
in the following theorem with the assumption that the mobile 
node does not have pause time. 

Theorem 1. Assume that S n = 0, Vro, and that Xq is at the 
origin. Then the spatial node distribution between Xq and X\ 
is characterized by the pdf f(r, 9) given by 



/M) = ^exp(-A7rr 2 ) 
7rr 



(10) 



Proof: See Appendix lAl ■ 
The physical interpretation of f(r, 9) is as follows. Let 
dA(r, 9) be a small area around the point (r, 9) given in polar 
coordinate. Then the probability P(dA(r, 9)) that the moving 
node resides in some measurable set A during the movement 
from Xq to Xi is approximately given by 



P{dA{r,6)) w f(r,e) ■ \dA(r,0)\, 



(11) 



where |cL4(r, 9) \ denotes the area of the set dA(r, 9). Also, as 
noted in the proof in Appendix lAl f(r, 9) can be regarded as 
the ratio of the expected proportion of transition time in the 
set dA{r,9) to |cL4(r,0)|. 

Now let us consider the case where the mobile node has 
random pause time S, characterized by the probability measure 
Pg(-). In this case, the spatial node distribution is given in the 
following theorem. 

Theorem 2. Assume that Xq is at the origin. Then the spatial 
node distribution between Xq and Xi is characterized by the 
pdff(r,9): 

f(r, 9) — p ■ exp(— A7rr 2 ) + (1 — p) ■ A exp(— A7rr 2 ), 
7rr 

(12) 



E[T] 



where p = 

y E[T} + E\S] 
tion time from Xq to Xi. 



is the expected proportion of transi- 



Proof: Note that the pause time S is independent of both 
X and V. Besides, the non-static probability p is the expected 
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proportion of the time that the node i is moving, i.e., 

El^tm E[T] 



a = 0.1 



P - 



lim ^fi 7T 



E[T] + E[S] ■ 



(13) 



The pdf f(r, 9) is the weighted superposition of two indepen- 
dent components as [27|: 



(14) 



where f(r, 9) is the spatial node distribution without pause 
time and is given in Theorem [TJ, and fx ^,9) denotes the 
spatial distribution of the waypoint X\ and is given in Lemma 
Q] Plugging f(r, 9) and fx ± (r, 9) into (TBI) , the expression ([T2l 
follows. ■ 

E. Model comparison 

In this section, we perform simulations to study the dif- 
ference between the proposed RWP mobility model and the 
classical one. For validation, we also compare them to a 
synthetic truncated Levy walk model [6|, which is constructed 
from real mobility trajectories. Thus, we indirectly compare 
the two RWP mobility models to real human walks. In 
the truncated Levy walk model, transition lengths have an 
inverse power-law distribution: Pl{1) ~ ttt^,0 < a < 2. 
The pause times also have an inverse power-law distribution: 
Ps(s)~-^,0</3<2. 

We simulate the movement of a mobile node under the 
three mobility models, respectively. The simulation area is a 
lOOOx 1000 gr idQ Since the main assumption in the proposed 
model is the Rayleigh distributed transition lengths, we first 
compare the statistics of transition lengths under the three 
mobility models. The results are shown in Fig. |2] where 
A is chosen such that the proposed RWP mobility model 
has the same mean transition length as that of Levy walk 
model. As shown in Fig. [2] transition lengths of the classical 
RWP mobility model are statistically much larger than those 
of Levy walk model. In contrast, transition lengths of the 
proposed RWP mobility model match those of Levy walk 
model better, especially in the low transition length regime. 
However, the proposed RWP mobility model lacks the heavy 
tail characteristic of the Levy walk model. This mismatch in 
the high transition length regime gradually diminishes as a 
increases. 

We next compare the statistics of direction switch rate under 
the 3 mobility models. The results are shown in Fig. [3] where A 
is chosen such that the proposed RWP mobility model has the 
same mean transition length as that of Levy walk model. The 
pause time distribution used in the 2 RWP mobility models are 
the same power-law distribution as that used for pause time 
distribution in the Levy walk model. As expected, the classical 
RWP mobility model has much lower direction switch rate in 
all the 3 subplots than the Levy walk model. In contrast, the 
difference in direction switch rate between the proposed RWP 
mobility model and Levy walk model is moderate when (3 
is large. When j3 is small, e.g., (3 — 0.1, the statistics of the 

'Note that, though the proposed RWP mobility model is defined on the 
entire plane for analytical tractability, in a simulation one has to deal with the 
boundary issue. Here we use reflection method. 
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Fig. 2. Simulated statistics of transition lengths: V = 1 and /3 = 1. 



direction switch rate of the proposed RWP mobility model are 
almost indistinguishable from those of the Levy walk model. 

The above comparison shows that the proposed RWP model 
is more flexible for mobility simulation in future cellular 
networks which have ever-smaller cells than the classical RWP 
model. In particular, some aspects of the proposed RWP model 
can be directly or indirectly fitted using traces, while the 
classical RWP model is fixed. Nevertheless, the comparison 
does not imply that the statistical difference between the 
proposed RWP and Levy walk model are non-significant. 
Also, we do not claim that the proposed RWP model can 
fully represent human walks. Also, human walks are not Levy 
walks even though the Levy walk model is constructed from 
real mobility trajectories. Indeed, human walks have complex 
temporal and spatial correlations and its nature has not yet 
been fully understood [7|. It is fair to say that none of the 
existing mobility models can fully represent human walks. 

Note that though many real mobility trajectories show that 
transition lengths of human walks seem to have an inverse 
power-law distribution, we do not choose it to model the 
transition lengths in our proposed RWP mobility model. The 
reason is that the mean transition length under inverse power- 
law distribution is infinite. This will cause analytical problems, 
e.g., the expected number of handovers in a movement period 
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Fig. 3. Simulated statistics of direction switch rate: V = 1 and a = 1. 



becomes infinite. Nevertheless, the analysis in this paper can 
be readily extended to any other transition length distribution 
which has finite mean. 

IV. Applications to Hexagonal Modeled Cellular 
Networks 

In this and the next section, we study the impact of 
mobility on important cellular network parameters, particularly 
handover rate and sojourn time, using the proposed RWP 
mobility model. We focus on hexagonal cellular networks in 
this section. 

A. Handover Rate 

We assume the typical mobile user is located at the origin. 
Then the expected number of handovers during one movement 
period can be computed as 



oo 



n=l , ' Cn 



(15) 



where P(dA(r, 9)) is the probability distribution of the way- 
point density fx x (r, 0) given in Lemma Q] and C n denotes 
the area covered by the n-th layer neighbouring cells. Now 
we formally define the handover rate. 



Definition 1. The handover rate is defined as the expected num- 
ber E[N] of handovers during one movement period divided 
by the expected period time. Mathematically, handover rate is 
given by H = E[N]/E[T P ]. 

Note that E[T p ] = E[T] + E[S] where E[T] has been 
given in Prop. Q]and E[S] can be determined from the pause 
time distribution. However, exact computation of N by (TT3T > 
is tedious. Thus, we propose the following approximation 
formula 

°° r-2-n r(2n+l)R 

E[N] app = y2n / f Xl (r,9)rdrd8, (16) 

n=1 Jo J{2n-\)R 

where R = \/Cpn with C being the hexagonal cell size. In 
other words, we approximate the 71-th neighbouring layer by a 
ring with inner radius (2n — 1)R and outer radius (2n + l)R. 
This approximation captures the essence of (fT5t and allows us 
to derive closed form results on E[N] and the corresponding 
lower and upper bounds. 

Proposition 3. Let d be the side length of the hexagonal 
cell and X the mobility parameter. The approximation of the 
expected number of handovers during one movement period is 
given by 



E[N] app = J^exp 



3\/3 



(2n + l) 2 Xd 2 



n=0 



and is bounded as E[N]% pp < E[N] app < E[N]% pp where 

e; y;; pp : . , -o(\ 



eVsxd 2 



6V3Xd 2 



1 - Q ^3V3Xd 2 



(17) 



(18) 



(19) 



Moreover, the difference AN app (Xd 2 ) between the upperbound 
and lower bound is a strictly increasing function of Xd 2 
and is within the range (0,1). In particular, AN app (Xd 2 ) — > 
as Xd 2 -> 00, and AN app (Xd 2 ) -> 1 as Xd 2 -> 0. 

Proof: See Appendix |B1 ■ 
Using Prop. [3] the approximate handover rate can be 
computed as iJ app = E } T , • E[N] a?p and is bounded as 



E\T„ 



where'lf^ = ^ • E[N] 
E[N]U As the size of the cells in cellular networks 



H.^ < i?app < H aw 



l H u 

app ? app 



becomes smaller, higher handover rates are expected. In this 
regard, it is interesting to examine the asymptotic property of 
handover rate as in Corollary |2] 

Corollary 2. Assume that any of the following asymptotic 
conditions holds: 1) d — > with fixed X; 2) X — > with fixed d; 
3) Xd 2 — > 0. Then the asymptotic approximate handover rate is 
given by 

I = I 

(20) 



app 



E[T] + E[S] V 6V3A 2d 



Though derivation by ( fl3T l is not tractable, the exact han- 
dover rate in hexagonal model can be obtained using a 
generalized solution of Buffon's needle [34| as in the following 
Proposition. 

Proposition 4. Let d be the side length of the hexagonal 
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cell and A the mobility parameter. The expected number of 
handovers E[N] during one movement period is given by 



E[N] = E[T] 



AV3E[V] 



37T d 

The handoverrate H = E[N]/E[T P ] is then given by 



H 



E[T] 4^/3 E[V] 



E[T] + E[S] 3vr d 



(21) 



(22) 



Proof: See Appendix [C] ■ 
We remark that the insights obtained by either approximate 
or exact approach are the same. Let us consider the simplified 
RWP mobility model where the mobile nodes do not have 
pause time and constant velocity, i.e., V = v, the asymptotic 
handover rate in Corollary |2] can be further simplified as 



7T V 



eVzd' 

and the exact handover rate H is given by 

3tt d : 



H 



(23) 



(24) 



which is consistent with the one given in [!35 ] . Note that the 
hexagonal cell size sh is given by S^fid 2 /2 in hexagonal 
tiling. So ( 1231 can be written as H Rpp ~ yT^-j==. Similarly, 



the exact method yields H 



Both results imply 



2 ^s H - 

that handover rate is inversely proportional to the square 
root of cell size sh- In other words, if we deploy more 
small cells and increase the BS density say by 4 times in 
current cellular network, we would expect the handover rate 
to be roughly increased by 2 times. Interestingly, the mobility 
parameter A does not play a role, while the velocity and the 
cell size affect the handover rate in a trade-off manner. Fig . [4] 
compares the number of handovers obtained by simulation to 
the counterparts evaluated by analytic formula dTTt and (l2ll 
respectively. It is shown that the exact analytic result closely 
matches the simulation while the approximation approach 
tends to underestimate the real number of handovers. 



B. Sojourn Time 

Sojourn time is defined as the expected duration that the 
mobile node stays within a particular serving cell. For brevity, 
we only consider the simplified RWP mobility model where 
the mobile nodes do not have pause time and constant velocity, 
i.e., V = v . Then sojourn time in any cell with coverage area 
C can be computed as 



S T = E[T]- / P(dA(r,0)) 



(25) 



where P(dA(r, 8)) is the probability distribution of the spatial 
node density f(r, 9) given in Theorem Q] For brevity, we only 
focus on the sojourn time in the cell where the connection is 
initiated during one movement period in the sequel. 

Assuming the mobile node is co-located with its currently 
associated BS at the origin, the sojourn time can be computed 



Simulation (d = 5) 
Simulation (d = 10) 
Simulation (d = 20) 
Exact Analysis (d = 5) 
Exact Analysis (d = 1 0) 
Exact Analysis (d = 20) 
Approximate Analysis (d = 5) 
Approximate Analysis (d = 10) 
Approximate Analysis (d = 20) 




4 5 

Mobility Parameter (X) 



Fig. 4. Handover in hexagonal model by analysis and simulation with velocity 
v = 1 and no pause time. 



as 



S T = E[T].^( 



d/2 r-y^d/2 



d/2 JO 



f(x,y)dydx 

— \/Zx + v3 d 



f(x,y)dydx), (26) 



where f(x,y) — exp(— A7r(a; 2 + y 2 ))/ \/x 2 + y 2 . However, 
no closed form result is available for (l26l . though it can 
be evaluated using numerical methods. We provide explicit 
formulas for the lower and upper bounds of (l26l in Prop. Q 
whose proof is omitted for brevity. 

Proposition 5. The sojoum time St can be bounded as < 

S T < S^ where 



S T *±E[T].(l-2QUI~n\d)) 



S^ = E[T] ■ (1 - 2Q(V27rAd)). 



(27) 
(28) 



Assuming that E[T] is finit^, St tends to since both 
the lower bound Sj- and upper bound Sj. tend to as d —> 0. 
This is an intuitive result. So cellular networks with small cells 
need to be equipped with fast algorithms. Otherwise, other 
handover strategies may be needed if the sojourn time in a 
cell is shorter than the time needed to complete the handover 
procedure. The analysis of the impact of mobility parameter 
A is more involved. To obtain insight, consider the constant 
velocity case which yields the following result. 

Proposition 6. Assume that V = v. Then as A —> 0, sojourn 
time St ~ otd/v, where a G 1) is a constant. 

Proof: See Appendix |Pl ■ 
From Prop. [6] we observe the interesting fact that sojoum 
time converges to ud/v as the mobility parameter A goes 
to zero. Then what matters is the velocity: Sojoum time is 

2 This is true for both the case with constant velocity and the case with 
uniformly distributed velocity on [v m - ln , v ma x] , as shown in Section IIII-BI 



s 



inversely proportional to the velocity. Also, we can express St 
as St ~ y^su ' \/^h wnere recall that sh is the hexagonal 
cell size. This shows that sojourn time is proportional to the 
square root of cell size, which is contrary to the asymptotic 
result for handover rate. 

V. Applications to Poisson-Voronoi Modeled 
Cellular Networks 

In this subsection, we apply the proposed RWP mobility 
model to analyze handover rate and sojourn time in cellular 
networks modeled by Poisson-Voronoi tessellation [1321 . Il36l . 
We first give a brief description on Poisson-Voronoi tessella- 
tion. Consider a locally finite set <f> — {cc;} of points Xi £ K 2 , 
referred to as nuclei. The Voronoi cell C Xi {4>) of point Xi with 
respect to 4> is defined as 



y - Xi || 2 < || y - xj \\ 2 ,yxj e (f>}. 



Let e x be the Dirac measure at x, i.e., for A £ 



,{A) = 1 



if a; £ A, and otherwise. Then the spatial point process 
<E> can be written as $ = 2~^ e ^i' an d trie Poisson-Voronoi 
tessellation is defined as follows 



Definition 2. For a spatial Poisson point process $ = J^. e Xi 
on M 2 , the union of the associated Voronoi cells, i.e., V( < &) = 
Lb e* ^i(^)' ls ca ^ e< ^ Poisson-Voronoi tessellation. 

In cellular networks modeled by a Poisson-Voronoi tessel- 
lation, the BSs are the nuclei distributed according to some 
PPP $ in M 2 . Besides, each BS Xi serves mobile users 
which are located within its Voronoi cell C Xi (<S>). The latter 
assumption is equivalent to the hypothesis of the nearest BS 
association strategy. In the sequel, we also assume that the PPP 
$ modeling the BS positions is homogeneous and its intensity 
is denoted by /i. 

A. Handover Rate 

Assume that the mobile node is located at the origin and let 
Xo and X\ be two successive waypoints. Conditioned on the 
position of Xi and a given realization of the Poisson-Voronoi 
tessellation, the number of handovers equals the number of 
intersections of the segment [X , Xi] and the boundary of the 
Poisson-Voronoi tessellation. Then we can obtain the expected 
number of handovers by averaging over the spatial distribution 
of Xi and the distribution of Poisson-Voronoi tessellation. 
This is the main idea used in proving Prop. [7j 

Proposition 7. Let /i be the intensity of the homogeneously 
PPP distributed BSs and A the mobility parameter. The expected 
number of handovers E[N] during one movement period is 

given by 

2 



E[N} = 



(29) 



The handoverrate H = E[N]/E[T p ] is then given by 

i 2 m 



H = 



E[T] + E[S] 7T 



(30) 



If we assume no pause time and constant velocity v, H in 
Prop. [7] can be simplified as 



(31) 



which is consistent with the one given in lF3~7l . In Poisson- 
Voronoi tessellation with the nuclei being homogeneous PPP 
$ of intensity fi, the expected value of the size sp of a typical 
Voronoi cell is given by sp = E{\C (&)\ 2 ] = 1/m EH, where 
C (<1>) is the typical cell and | C ( ) 1 2 denotes the area of 
C (<1>). So if we assume no pause time and constant velocity 

4 

v , H in Prop. |7jcan be simplified as H = —v/yfsp, which 

7T 

implies that the handover rate is inversely proportional to the 
square root of the cell size. This is consistent with the results 
in the hexagonal model. 

Fig. [5] illustrates that the analytical result ((29) matches the 
simulation result quite well. Also, we compare the handover 
rate of Poisson-Voronoi model, exact and approximate han- 
dover rate of hexagonal model in Fig. [6] as a function BS 
intensity. They all indicate that handover rate grows linearly 
with the square root of the BS's intensity ^fjl. We further eval- 
uate the three types of analytic results on handover, i.e., Prop. 
[3 Prop. |4] and [7j by simulating the proposed RWP mobility 
model using the real- world data of macro-BS deployment in a 
cellular network, provided by a major service provider. Recall 
that we assume each BS serves the mobile users located within 
its Voronoi cell. A handover occurs when the mobile crosses 
the cell boundary. Thus, only the BS location data are relevant 
in the simulation. There are 400 BSs which are distributed in 
a relatively flat urban area. These 400 BSs roughly occupy a 
105 x 90 km area. We normalize the network size to be 1 x 1 . So 
/i = 400 in the Poisson-Voronoi model, while the side length d 
in hexagonal model is determined through 3\/3d 2 /2 — 1/400. 
The results are shown in Fig. [7J It can be seen that the 
Poisson-Voronoi model is about as accurate as the hexagonal 
model in predicting the number of handovers. Meanwhile, 
the approximate analytic result underestimates the number of 
handovers. 



B. Sojourn Time 

Unlike the hexagonal model, (f25T > is not sufficient for 
computing the sojourn time in Poisson-Voronoi tessellation 
modeled cellular networks. Indeed, assuming that the mobile 
user is located at the origin, the sojourn time in the cell 
involves an additional source of randomness - the Poisson- 
Voronoi tessellation. So even averaging over the spatial node 
distribution sojourn time is still a random variable. In the 
sequel, we aim to characterize the probability distribution of 
this sojourn time. To this end, we first introduce the concept 
of contact distribution. Consider a random closed set Z and a 
convex compact test set B containing the origin o. Then the 
contact distribution is defined as |29| 



H B (r) = P(ZnrB ^ 0|o ^ Z),r > 0. 



(32) 



Proof: See Appendix [E] 



In this paper, we are particularly interested in a special case: 
the linear contact distribution function Hi(r) with the test set 
B being a segment of unit length /. Since Poisson-Voronoi 
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Fig. 5. Handover in Poisson-Voronoi model by analysis and simulation with 
velocity v = 1 and no pause time. The number of handovers is proportional 
to the square root of the BS's intensity /i, and is inversely proportional to 
mobility parameter A. 



Fig. 7. Evaluation of handover by an actual macro-BS deployment with 
velocity v = 1, no pause time and normalized BS's intensity fj, = 400. 
Poisson-Voronoi model is about as accurate in terms of handover evaluation 
as hexagonal model. 



* Poisson-Voronoi handover rate 
□ Hexagonal exact handover rate 

Hexagonal approximate handover rate 




0.05 1 1 1 1 1 1 1 1 ' 

0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 

Intensity of Base Station (p.) 

Fig. 6. Comparison of handover rate in Poisson-Voronoi model, exact 
and approximate handover rate in hexagonal model with velocity v = 1, 
normalized BS's intensity /i = 1, mobility parameter A = 1 and no pause 
time. They all imply that the handover rate grows linearly with ^//i. 



tessellation is isotropic, the orientation of this segment I is not 
important. For Poisson-Voronoi tessellation modeled cellular 
network, the random closed set Z of interest is the union of 
all cell boundaries U Xie $9C Xi (<&), where dC Xi (<£>) denotes the 
boundary of the cell C Xi ($). 

With a slight abuse of notation, let C ($) denote the Voronoi 
cell containing the origin o. We can further simplify Hi (r) as 
follows: 

P(Z n rl -- 



Hi(r) = 1 



= 1 - P(rl C C„(#)). (33) 



1 - P(o e Z) 

The last equality follows because 1) the origin o is contained 
in the interior of C ($) almost surely and thus P(o £ Z) = 
where Z = U Xie <s,dC Xi (3>), and 2) the event {Z n rl — 0} is 



equivalent to that {rl C C„($)}. Now we are in a position to 
characterize the sojourn time St- 

Proposition 8. Let p be the intensity of the homogeneously 
PPP distributed BSs and A the mobility parameter. The pdf of 
sojourn time St is given by 



fsAt) 



2y^XE[T] V 2 



E[T\ 



W)hi (x) 



(34) 



where E[T] is the expected transition time, x 

r) 

sin 2 a sin (3 



J=Q(-D(I(1 - ^)) andh^r) is given by 



h[(r) = 4717/ 



7T P7T — 0. 



MP) 



jo sin (a + /3) 

exp(-pV 2 {r, a, /3))dj3da, r > 0, (35) 



V 2 {r,a,f3) = 7rp 2 (a (a) + ai(a)) + 7r(r 2 + p 2 - 
>a)(a (/3) + aM) with p = ™^ >ao (6) = 1 - &, 



where Vi(r, a, f3) 

2rpcosa)(a (p) - l „,u, H - sin(Q , +j3) 

and ai {0) = ^M,wdb (p) = ( ^ /3)co ^ +sin/3 . 
Proof: By Theorem [U we have 

r-2-K r-r 

Sr(r) = E[T] / / cxp(— \irx 2 )xdxda 

Jo Jo nx 

= E[T] • (1 - 2Q(v / 2rrAr)). 

Clearly, F ST (t) = if t < 0, and Fg T (t) = 1 if t > E[T]. If 
0<t< E[T], 

F ST (t) = P{S T <t)= P{E[T] • (1 - 2Q(v / 2 T rAr)) < t) 



= P[ r 

= H 



E[T] 



)) 



'(Ki-irfrr)) 



hi(r)dr. 
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Fig. 8. Plot of pdf of sojourn time in Poisson-Voronoi model with mobility 
parameter A = 0.01 and velocity v = 1. 



Taking the derivative with respect to t for both sides yields 
( |341 l. The closed form expression of hi(r) in ( f33T > has been 
given in |38|. This completes the proof. ■ 

Using the previous results for expected transition time, the 
pdf of the sojourn time under different velocity distributions 
can be derived from Prop. [8] Note that the distribution of 
the sojourn time is instrumental in studying the ping-pong 
behavior for mobility enhancement in cellular networks |H]. 
Though a closed form expression is not available, the pdf of 
the sojourn time given in Proposition |8]only involves a double 
integral, which is tractable for efficient numerical evaluation. 
We plot the pdf of sojourn time in Fig. [8] for illustration. As 
expected, the smaller /i is, the more likely the mobile node 
stays longer within the cell. This is intuitive since a smaller 
/i implies larger cell sizes on average. As a result, it is less 
likely that the mobile node would move out of the cell. 

We also compare the analytic result about sojourn time in 
the Poisson-Voronoi model to its deterministic counterpart in 
the hexagonal model in Fig. [9] It is shown that the analytic 
result in the Poisson-Voronoi model is more conservative and 
yields smaller mean sojourn times than its counterpart in the 
hexagonal model. Besides, we can see that the upper and lower 
bounds of the sojourn time in hexagonal model are pretty tight. 

VI. Conclusions and Future Work 

In this paper, we study the critical mobility issue in cellular 
networks. To this end, we first propose a tractable RWP 
mobility model defined on the entire plane. The various 
properties of the mobility model are carefully studied and 
simple analytical expressions are obtained. Then we utilize this 
tractable mobility model to analyze the handover rate and so- 
journ time in cellular networks. The analysis is carried out for 
cellular networks under both hexagonal and Poisson-Voronoi 
models. We derive closed form expressions and/or bounds for 
the performance metrics in question. These analytical results 
are instrumental for mobility management in cellular networks. 



Fig. 9. Comparison of sojourn time in hexagonal tiling and Poisson-Voronoi 
tessellation with velocity u = 1. Hexagonal model yields larger sojourn time 
than that of Poisson-Voronoi model. Besides, the upper and lower bounds for 
hexagonal model are tight. 



Note that the proposed RWP mobility model represents the 
real movement of mobile nodes in a simplified manner. Thus, 
it does not capture some other mobility characteristics such 
as temporal and spatial dependency of the mobility pattern. It 
is desirable to extend the current model further to incorporate 
these extra mobility characteristics while maintaining a certain 
degree of tractability. It would be rather interesting (and 
challenging) to characterize the intuitive tradeoff between the 
complexity and tractability of the mobility models. Besides, 
it is also interesting to evaluate the performance of the many 
wireless protocols by applying the proposed tractable model 
in theory and/or simulation. 
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Appendix 

A. Proof of Theorem Q] 

We first derive the pdf of the random waypoint X\ = 
(Rl, Oi) as follows: 

P(Ri <r + Ar) - P{Ri < r) 



f Xl (r,0) = lim 



lim 

Ar-tO 

lim 

Ar->0 



exp(— Xirr 2 ) — cxp(— \ir{r + Ar) 2 ) 
2nrAr 

27rA(r + Ar) exp(-A7r(r + Ar) 2 ) 
2nr 



= Aexp(— Xiirr 2 ). 
This result is summarized in Lemma[TJfor ease of reference. 
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Lemma 1. Given Xo is at the origin, the pdf of the random 
waypointXi = ©i) is given by 

fx 1 {r,e)^Xexp(~Xirr 2 ). (36) 



We next derive the pdf of the spatial node distribution. The 
main technique of the following proof is inspired by [26 1, 
which is also adopted in l25ll . Consider a small set dA located 
at (r, 6). Let L\ denote the vector X\ — Xo and \L\ \ = L\. 
Then the pdf of the spatial node distribution can be interpreted 
as the ratio of the expected proportion of transition time in the 
set dA to the area |cL4|, i.e., 

Ej^ndAyv] _ E[\L 1 ndA\] 
J[ ' ' E[Li/V]|<L4| ^[£i]|<L4| ' K 

where the second equality follows from the independence of 
V and the waypoints. Note that 

E[\L X n dA\] _ C f Xl {x, 6)xdxd9 ■ Al 



\dA\ r ■ d9 ■ Al 

J r A exp(— \irx 2 )xdx 1 



= o — cx p(-Atit ), 
2irr 



(38) 



where Al denotes the length of the small intersection if L\ 
intersects dA, and we apply Lemma [TJ in the second equality 
in (1381 . The first equality in (1381 can be explained through Fig. 
[Tolas follows. The intersection of \L\ H dA\ is Al if X\ is in 
the shaded area and otherwise. Noting that the probability of 
the event that X\ is in the shaded area is J r °° fx t (x, 8)xdxd9, 
we have 

/>oo 

E[\L X n dA\] = Al ■ / f Xl (x,6)xdxde. (39) 

J r 

This and \dA\ = r • dO • Al give the desired equality. So 

f(r,0) = eX ^f ] ) = VXeM - X7tr2 \ (40) 
2i:rhj[Li\ nr 

where we use the result that = 1/2 -\/A in the last 

equality. This completes the proof. 

B. Proof of Proposition \3\ 

From Lemma [JJ /x[(r,ff) = A exp(— Airr 2 ). We compute 
E[N]spp as follows. 

E[N]. dpp = J2n / Xexp(~\Trr 2 )rdrdd 

n=1 JO J(2n-1)R 

oo 

= n (exp(-7r(2n - l) 2 i? 2 A) - exp(-7r(2n + l) 2 i? 2 A)) 

n=l 

oo 

= Y exp(-7r(2n + l) 2 i? 2 A) 



n=0 



[C tf27d 
Substituting R = \l — = — j=, we obtain 



7T V27T 



E[N], pp = exp 



3V3 



(2n + l) 2 d 2 X 



Then the lower bound A^„ mv is derived as follows. 

tip pi UA 



E[Nl 



> 



exp 



3^3 



(2x + l) 2 d 2 \ dx 



eVzxd 2 



Q JzVsXd 2 4 -E[iV]j 



Japp- 



The upper bound i?[iV]^ p can be derived in a similar manner. 

For the remaining proof, denote t = \J 3\/3Ad 2 . Note that 
t > 0. Then the difference AiV app = E[N]g p - E[N]^ pp can 
be compactly written as v/f j(l — 2Q(i)) which we denote 
by We claim <?(£) is strictly decreasing when t > 0. To 

see this, define h(t) = y^t exp(— -y) + 2Q(t) — 1 whose 

derivative is "— \J~^;t 2 exp(— —)" which is strictly negative 
when £ > 0, implying is strictly decreasing when t > 0. 
Meanwhile, h(0) = 0. So h(t) < when t > 0. Now it is 

clear that the derivative of given by c/ 1 ^ (t) — \ ^- is 

y 2 t 2 

strictly negative when i > 0. Thus, is strictly decreasing 
when t > 0. Besides, 

£2 

lim */--(! - 2Q(t)) = limexpf ) = 1, 



t-vo 
lim 

f— f oo 



~(1 -2Q(t)) = lim exp(- 



7T 1 



2 t 



(l-2Q(i)) = 0. 



The desired results follow by further observing t is a strictly 
increasing continuous function of Ad 2 . 

C. Proof of Proposition 

We use (generalized) argument for Buffon's needle problem 
in this proof. Consider the typical node located in area A of 
size \A\ and cell boundaries Ba of length \Ba\- Then the prob- 
ability that this node crosses the small boundary Ab within a 
small time interval At is Ap = -^-Ab- vAtE[\ sin 6|], where 
is uniformly distributed on [0, 2n] (following from our 
mobility model) and thus E[\ sin0|] = -. So the probability 
that this node crosses B A is p(A) = Ap ■ ^ = • vAt. 

Thus, conditioned on moving, the handover rate 



H* 



lim 

\A\-yo 



At 
9d 



—v ■ lim 

4V3^ 



\A\ 



tt 9V3d 2 /2 3vr d' 
Correspondingly, E[N] = E[H*\ ■ E[T] = ^M^3 . E[T] 



and H 



E[T] 4^3 E[V] 
E[T] + E[S] "3^ d~ ' 



3tt 



D. Proof of Proposition [6] 

Applying the general expression (l26l for St yields 



St = — f{£xd,£ y d)- 

TTIS 

d/2 r V5d/2 



Jo 



dydx 



d/2 JO 



-V3x+V3d 



dydx, 
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Fig. 10. A geometric illustration for the proof of Prop. [T] 



which equals 



3v / 3exp(-Avr(c:2d 2 + Qd 2 )) d 



2vr 



(41) 



where (£ x ,£y) are a pair of constants in the region A := 
{(x,y) : x G [0,1], y G [0, ^] if a; G [0, ±) and y G 
[0, \/3(l — £ x )] if x G an d in the last equality we 

use E[T] — - ^= derived for constant velocity case in 
Section [nil an d mean value theorem for integrals. Then, 



liniA^o S T 
constant a 
® 



3V3 



Next we give bounds on the 



2 ^ v^I+?I 

— 3Vji 1 Using the bounds given in Prop. 



2tt 



lim 5-^ = lim 

A->-0 A^O 



1 - 2Q (^f^Ad) 



2i/-\/A 



\/3d . 3 2 . \/3d 
lim exp — ixd A = . 



(42) 



Similarly, lim^o S 1 ^ = djv. Since these bounds are strict, 
we conclude that a G (-^,1). 

£. Proof of Proposition \7\ 

This proposition can be proved by following the proof of 
Prop. [4] Nevertheless, we provide an alternative proof here. 
To this end, we first introduce some terminologies for ease 
of exposition. Consider the set x^ consisting of arbitrary 
k distinct points in $ on K 2 . Without loss of generality, 
we assume x^ = {xi, Xf-}. Then if the intersection 
T(x^\<S>) = n^ =1 C Xl ($) ^ 0, then J"(a;( fc )|$) is called a 
(3 — fc)-facet. Now let be the set of all configurations 
x (3-k) gj <jj That is, the intersection of the Voronoi cells 
associated with any configuration G is a fc-facet. 

For each configuration a;( 3 ~ fe ) = {x%, x^-h} G $fe, we 
associate a point z x (3-k) ($) called "centroid" such that for all 



y G K 2 , z(x 



(3-fe) 



y|$ + y) = y + z(cc (3 ~ fe) |$). Note that 



there are several degrees of freedom in choosing the centroids. 
We refer to ll36l for more details. 

Consider the intersection of the Poisson- Voronoi tessellation 
with a fixed line L. Without loss of generality, we assume 
that L contains the origin. Then the nonempty sectional cells 
C*,($) = C Xs ($) n L satisfy that (i) L = U XtE ^C x , ($), 
and (ii) ri(C Xi {$)) H ri(C Xj ($)) = 0,Vi ^ j where ri(-) 
denotes the relative interior [36 1 . So the sectional cells C Xi (<£>) 



constitute a tessellation of L, denoted by Vl($). Now consider 
the intersection T L (xW\$) — F(x^\<&) n L which can be 
either empty or a singleton. Clearly, 7-"i(cc^|$) is the 0-facet 
of tessellation Vl(<I>) if Fl{x^ 1$) is nonempty. Now let 
*J = {x^ = {xi,x 2 } G $1 : T L (xW\$) ± 0} which 
parametrizes the 0-facets of Vl($). Then the intensity of the 
0-facets of Vl(3>) is well defined as 



Mo 



E[E x ^^x{z L (x^meB L }} 
\Bl\i 



(43) 



where x{'} lS indicator function taking value 1 if the event 
in its argument is true and otherwise, zl(x^\^>) is the 
centroid of the 0-facets of Vl($), Bl Q L is any arbitrary 
Borel set and \Bl\i is the volume of Bl with respect to K. 

Without loss of generality, assume that Xo is at the origin. 
Let L := L(Xq,Xi) be the line containing Xo and X\, 
and Bl the interval [Xo, Xi]. Then conditioned on the next 
waypoint Xi = (r, 8), the expected number of handovers can 
be computed as 

E[N\X 1 = (r,e)]=E[ x{xeB L }} 
xeu a . i e$dc a . i (*)nL 

= E[J2 x{zl(xW\$) G B l }} (44) 
4r^/J2 



Mo l-^o — Xi|i — 



(45) 



where dC Xi (&) in (l44b denotes the boundary of the cell 
C Xi ($). The equality in ( f44b follows by choosing the centroids 
zl(x^ as follows: for any x^ G choose the singleton 
.F(a: (2) |<p) n L as the centroid ;z L (a: (2) |$) of the 0-facets 
of Vl( < 1 ) ). The equality in ([45]) follows from (|43j. The last 
equality follows since [i$ — A^Jl/tt for Poisson- Voronoi 
tessellation with intensity /i [36|. Then handover rate can be 
computed as follows: 

E[N] E[E[N\X 1 = (r,0)]] 



H 



E[T] E[T] 

j />2tt poo 

—jrj^ j o S[AT|X 1 = (r,e)]/x 1 (r ! %drdrJ 
1 f 2 * r°4ry/I, 



E[T] 

i 4r(|) 

E[T] n i 
1 2 //I 



-Aexp(— \itr )rdrd9 



(46) 



(47) 



S[T] 7T V A' 

where in ( l46b we use Lemma [T] and in ( f47T > we apply the 
result that J °° r Q_1 exp(— ^nr^) dr = for a, /3, 7 > 0. 

Plugging -E[T] yields the desired results. 
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